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Abstract 



The non-equilibrium dynamics of a Yukawa theory with N fermions coupled to a scalar field is studied 
in the large N limit with the goal of comparing the dynamics predicted from the renormalization group 
improved effective potential to that obtained including the fermionic backreaction. The effective potential 
is of the Coleman- Weinberg type. Its renormalization group improvement is unbounded from below and 
features a Landau pole. When viewed self-consistently, the initial time singularity does not arise. The 
i-S^ ' different regimes of the dynamics of the fully renormalized theory are studied both analytically and 

numerically. Despite the existence of a Landau pole in the model, the dynamics of the mean field is 
^ • smooth as it passes the location of the pole. This is a consequence of a remarkable cancellation between 

^ I the effective potential and the dynamical chiral condensate. The asymptotic evolution is effectively 

■ described by a quartic upright effective potential. In all regimes, profuse particle production results in 
the formation of a dense fermionic plasma with occupation numbers nearly saturated up to a scale of the 
order of the mean field. This can be interpreted as a chemical potential. We discuss the implications of 

■ these results for cosmological preheating. 

^' 

1 Introduction 

In the Standard model or its extensions, the Yukawa couplings of Fermions to scalars (Higgs) play a funda- 
mental role. Not only do such couplings determine the masses of the fermionic degrees of freedom, but in 
turn, it is through these couplings that the fermionic sector influences the dynamics of the scalar fields. If 
these Yukawa couplings are large enough, they can lead to negative contributions to the beta functions of the 
running scalar self-couplings and so to destabilizing the vacuum by large negative radiative corrections to the 
scalar eS'ective potential Q]. This is the Coleman- Weinberg mechanism of symmetry breaking by radiative 
corrections 0. 
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While such a scenario has been ruled out within the Standard model due to an unacceptably low value 
of the Higgs and the top quark masses, large negative radiative corrections to the effective potential from 
large Yukawa couplings could still be relevant in extensions of the Standard model with more complicated 
Higgs- Yukawa sectors]^. Coleman Weinberg phase transitions in extended Higgs models and their potential 
cosmological implications have been studied by Sher||^ who analyzed the effective potential of this theory. 
While this study extracted bounds on the parameters of extended Higgs sectors from vacuum stability and 
thermodynamic considerations, these results are based on an equilibrium description based on the effective 
potential, a purely static quantity. 

Detailed studies reveal that the information extracted from a static effective potential is restricted to 
situations very close to equilibrium and that a deeper understanding of dynamical processes requires a non- 
equilibrium treatment. The necessity for a non-equilibrium description of quantum field theory has become 
clear in cosmology where inflationary phase transitions require a fully dynamical description [Q, in heavy- 
ion collisions where a transient quark-gluon plasma may be formed ^ and domain formation in phase 
transitions 10, || which may have consequences in cosmology as well as in heavy ion collisions. 

In each of these fields, non-equilibrium effects can give rise to new phenomena which can differ from 
equilibrium behavior in dramatic ways that cannot be captured by an effective potential description. 

From the point of view of cosmology and inflationary cosmology in particular, non-equilibrium effects 
associated with particle production via resonances and or instabilities in bosonic field theories have taken 
center stage. This is evident in the theory of preheating|^ as well as in the classicalization of fluctuations 
during inflation |[fC||, where spinodally unstable dynamical fluctuations about an homogeneous condensate 
modify the long wavelength behavior of the theory. 

While the non-equilibrium dynamics of fermionic fields has been studied recently by several authors 
p] , |T2|, Q and while there has been some remarkable progress in lattice simulations of fermionic dynamics 
in low dimensional gauge theories the dynamical behavior of fermions has not received the same level 
of attention as the bosonic case. This is mostly due to the argument that Fermi-Dirac statistics preclude 
parametric amplification of occupation numbers with non-perturbative particle production. 

In order to have access to non-perturbative dynamics, in this work we consider theories containing N 
fermions with Yukawa couplings to one scalar field $ in leading order in the large N limit. The large N 
limit leads to a consistent, non-perturbative approximation scheme that can be systematically improved. 
This approximation has been applied to bosonic theories and has been studied both analytically to leading 
order [l^ and more recently numerically in real time including corrections beyond leading order in the 
large iV|l7|, [l^. These studies reveal a wealth of new phenomena not accessible via perturbative methods. 

In the large N limit we study here, the fermions serve to suppress scalar field fluctuations to leading 
order in 1/N. We find that even at leading order a new and important aspect of field theory comes into 
play: the wave function renormalization. It should be noted that this does not occur at the same order in 
purely scalar theories. Our study reveals that this new ingredient is responsible for dramatic new dynamical 
phenomena. 

We obtain the equations of motion of the mean field, or expectation value of the scalar field, taking into 
account the non-equilibrium back reaction of the fermionic modes to leading order in This analysis 

goes beyond the effective potential approximation!^ and affords us a window to understanding the non- 
perturbative dynamics of this coupled system. 

We can summarize our main results as follows: 

• The fully renormalized theory displays the Coleman- Weinberg mechanism of dimensional transmutation 
through radiative corrections. It has an effective potential that is unbounded from below at large values 
of the mean field with two symmetric global maxima and a local minimum at the origin. It also exhibits 
a Landau pole at an energy scale which is non-perturbative in the Yukawa coupling. The presence 
of the Landau pole distinguishes two distinct regimes to be studied depending on the relationship 
between the cutoff A, which is required for the numerical analysis of the theory, and the position 

of the Landau pole. 

• Suppose that we take A < E^. Then, if the initial value of the mean field is between the origin and the 
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global maxima, the mean field oscillates about the origin, and the fcrmionic quantum fluctuations grow 
as a result of particle production in a preferred band of wave vectors. This is akin to what happens in 
the bosonic case but here the production must saturate due to Pauli blocking. The width of the band 
of wave vectors is determined by the initial amplitude of the mean field and the mass of the scalar field. 
If, however, the initial value of the mean field is larger than the maxima, its amplitude runs away to 
the cutoff scale at which point the evolution must be stopped since the theory reaches the edge of its 
domain of validity. 

• Our most noteworthy results are those for which A >> El and the initial value of the mean field is 
taken to be larger than the position of the maxima of the effective potential. In this case, an analysis 
of the dynamics based solely on the static effective potential would lead to the conclusion that the 
time evolution of the mean field would lead to a divergence or discontinuity in the time derivatives 
when the amplitude reaches the value of the Landau pole. However a detailed analysis of the full 
dynamics, including the fermionic fluctuations and their backreaction on to the mean field, reveals 
that the evolution of the mean field is sm,ooth. When the amplitude of the mean field becomes larger 
than the Landau pole, the dynamics becomes oscillatory and asymptotically reaches a fixed point 
described by a simple quartic, upright effective potential with a quartic coupling of order one and with 
the mean field oscillating with a large amplitude around the origin. This novel dynamical behavior 
arises from a remarkable cancellation between the fermionic fiuctuations and the contribution from the 
instantaneous effective potential that leads to smooth dynamics through the Landau pole. When the 
initial value of the mean field is larger than the maxima of the effective; potc^ntial but much smaller 
than the Landau pole, the ensuing non-equilibrium evolution leads to fermion production in a band 
of wave-vectors up to the scale of the Landau pole. These modes become populated with almost 
Pauli blocking saturation at large times and describe a very dense medium. We study this behavior 
numerically and confirm that this phenomenon occurs for a wide range of parameters. We are led to 
conjecture that the theory is, in fact, sensibly behaved beyond the Landau pole when studied both 
dynamically and non-perturhatively, at least at the mean-field level. While at this point this is merely 
a conjecture, these phenomena may have some interesting phenomenological consequences. 

• A consistent analysis of the renormalization during the dynamical evolution reveals that the wave func- 
tion renormalization builds up in time over a time scale of 0(1/A). The fully renormalized equations 
of motion display a renormalized coupling at a scale determined by the amplitude of the mean field, 
and which therefore depends parametrically on time. This is a consequence of the "running" of the 
coupling constant with scale, which in the dynamical evolution translates to a "running" with time. 

The article is organized as follows: in section 2 we obtain the renormalization group improved effective 

potential, discuss its features including the presence of the Landau pole and the potential singularities that 
would occur in an analysis based solely on the effective potential. In section 3 we obtain the equations of 
motion to leading order in the large N limit. In section 4 we address the renormalization of the equations 
of motion and the energy density. We discuss and resolve the issue of potential initial time singularities; in 
particular, we highlight the fact that the wave function renormalization builds up on time scales determined 
by the cutoff. In this section we also establish contact between the non-equilibrium equations of motion and 
the renormalization group improved effective potential, emphasizing the emergence of smooth dynamics as 
the mean-field approaches and passes the Landau pole. In section 5 we provide a detailed and comprehensive 
numerical study of the dynamics in several cases in a wide range of parameters. In section 6 we provide an 
exact proof of the lack of unstable bands for fermionic mode functions in the background of a scalar field 
that oscillates, by showing that the Floquet indices are purely real. This is the underlying reason for Pauli 
blocking at the level of mode functions. We also offer a perturbative analysis of this important phenomenon, 
which provides the reason for the existence of a preferred band of wave vectors for the produced fermions. 
We summarize our findings and offer some conjectures for potential implications of our results in cosmology 
as well as for the phenomenology of theories with extended Higgs- Yukawa sectors containing heavy fermions 
(and hence large Yukawa couplings) that could feature a Landau pole in an energy range of phenomenological 
interest. 
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2 Static aspects: the effective potential and its RG improvement 



Before studying the dynamical aspects of the Yukawa theory in the large N limit it is illuminating to 
understand the static aspects via the effective potential and its renormalization group improvement. 
The Yukawa model under consideration is described by the Lagrangian density 



2^^^ 2 4!iV 

N 

- Htd 

V-., (2.1) 



where the subscript B denotes the bare fields, anticipating the need for renormalization. The factors of 
N in the coupling constants are explicitly displayed so that both the quartic self coupling and the Yukawa 
coupling are of 0(1) in the large N limit. 

The effective potential is defined as the expectation value of the Hamiltonian in a state that minimizes 
the energy subject to the constraint that the field has a space-time independent expectation value p9| 

(^\H\^) 

ye//($)-min^^^ (2.2) 
with fl the spatial volume and |$) a (coherent) state for which 

(|.|$||.)=$. (2.3) 

To leading order in the large N the effective potential is obtained by replacing $(x) ^ $ in the Hamilto- 
nian and neglecting the scalar field fluctuations, since the energy will be dominated by the N fermion fields. 
This is the mean field approximation, which becomes exact in the large N limit. 

Hence 



Hm = n (\ml $2 + $4 V ^ . V + M] v., (2.4) 

M = (2.5) 

The fermionic contribution to the Hamiltonian is simply that of N Dirac fermions of mass M, and can be 
diagonalized in terms of creation and annihilation operators for particles and antiparticles with dispersion 
relation lu^ = \/W^+M^ . The state that minimizes the expectation value of the normal ordered Hamiltonian 
is the vacuum state for particles and antiparticles and corresponds to the Dirac sea completely filled (with 
two spin states per wave vector). It is convenient to rescale the expectation value of the scalar field to exhibit 
the large N limit more clearly: 



in terms of which 



^ = y/N 5b (2.6) 



Veim = ^ ^ ^-^l + ^4 - 2 / ^ V'fc^TlP \ . (2.7) 



t2 



The integral in (2.7) is the (negative) contribution from the Dirac sea, and is calculated with an ultraviolet 
cutoff A. A straightforward calculation, subtracting the "zero point energy" cx A^ and neglecting terms that 
vanish in the limit A — > oo, leads to 
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(2.8) 



where ji is an arbitrary renormalization scale. The renormaUzation of the mass and quartic scalar coupling 
A can be gleaned directly from the form of the effective potential above. However, to make contact with the 
dynamics in which the equations of motion are obtained from the non-equilibrium effective action, we also 
need the field or wave function renormalization. The wave function renormalization cannot be extracted 
from the effective potential since it is associated with gradient terms in the effective action. To leading order 
in it can be obtained from the one-fermion loop self energy shown in figure (|l|) below, where the scalar 
field lines have non-vanishing external momentum. 



Figure 1: Scalar self energy in leading order in large N . External dashed lines correspond to the scalar field, 
internal solid lines are fermions. 



This calculation leads to 



ZM=(l + ^ln 



2A' 



(2.9) 



The renormalization of field, mass and couplings is now achieved by introducing 
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in terms of which the renormalized effective potential is given by 



6i / A, 



Vr <5i?e4 



(2.10) 

(2.11) 

(2.12) 
(2.13) 

(2.14) 



We will consider the case in which the renormalized mass of the scalar field vanishes, since this case 
will highlight the important feature of dimensional transmutation both at the level of the static effective 
potential as well as the dynamics. Therefore in what follows we set ttib. = 0. 
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The equation of motion, which wiU be the focus of next section, requires V'{Sii) — dV{Sij) /OSr. This is 
given by 



where we introduced the renormaUzed effective ferniion mass 
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In terms of the scale M we find 



= 0. 



(2.15) 



Mr = yR^R. (2.16) 
We see that the effective potential features an extremum at Mr = M with M determined by 



(2.17) 



VeffiSR) 
VeffiSR.) 
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87r2 
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M 



(2.18) 
(2.19) 



Thus we see that the term Aij/3! inside the parenthesis in ( ^.14 ) above has been traded for a new 
scale M at which the effective potential features a maximum. This is the manifestation of dimensional 
transmutation . 

Figure d) below displays Veff/M^ vs. X = The 14// is given by (|t|). 



IL 




-1.5 -1.0 -0.5 0.0 0.5 1.0 



Figure 2: F{x) = -X^ Indxle"^) vs. x = ^ 



While there are alternative calculations of the effective potential, the Hamiltonian formulation highlights 
many important aspects that will be relevant to the discussion in the next sections. In particular, it makes 
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clear that the effective potential is the expectation value of the Hamiltonian in a "vacuum" state in which the 
scalar field attains an expectation value. Furthermore, this Hamiltonian interpretation immediately provides 
the physical reason for the effective potential being unbounded below in this approximation: it is completely 
determined by the negative energy Dirac sea. For larger amplitudes of the expectation value, the effective 
fermionic mass is larger, thus the negative energy of a free fermion mode of momentum k in the Dirac 
sea decreases further. Thus for larger amplitudes of the expectation value, the energy stored in the Dirac 
sea becomes more negative. This particular point should be borne in mind when we study the dynamics of 
the mean field below, since we will find that the evolution of the scalar field "feeds off" the negative energy 
Dirac sea. 



2.1 RG improvement: 

The emergence of the scale M is a consequence of the renormalization scale fi introduced above. A change 
in this scale is compensated for by a change in the couplings. The effective action 



(Tx 



(2.20) 



is invariant under a change of the renormalization scale, and consequently under a change in the scale M . 
This invariance leads to a renormalization group equation for the effective action 

Since our main focus is to study dynamical behavior, which involves Viff, we now use the renormalization 
group to improve the derivative of the effective potential. 

While in principle we can study the full solution of the renormalization group equation as in 1), the large 
N approximation simplifies the task. In this limit we need only keep the one loop fermion contribution to 
renormalization. Therefore, after trading the quartic self-coupling for the scale M via dimensional transmu- 
tation, the effective potential (and its derivatives) are only functions of the Yukawa coupling. Furthermore, 
from the renormalization conditions (2.1C,2.13) the product i/rSj^ is a renormalization group invariant, i.e, 
is constant under a change of scales. With the purpose of comparing with the dynamics to be studied in the 
next section, it is convenient to introduce the coupling 



9r(p) 



2^2 



and to RG improve the product 



yR{^^)viff{SR{^i)) = -gR{^i)Ml\l, 



Mr 

M 



(2.21) 



(2.22) 



The reason for studying this product is based on the idea that the effective equation of motion of the 
scalar field is loosely of the form 



5R{t) + Viff{5R{t))^Q. (2.23) 

While 5 is not invariant under a change of scale (i.e, under a renormalization group transformation) in 
the large A'^ limit the product Mr^ = UrSr is a renormalization group invariant. Thus one is led to consider 
the product yRV^jrj:. 



However, from the renormalization condition (2.12) of the Yukawa coupling and the wave-function renor- 
malization constant in the large N limit given by (2.9) we find 



2 2 , 

— - - — + In 

9r(P-) 9b 



2A 



which leads to the renormalization group running of this coupling 

2 2 



5fl(M') .9i?.(A*) 



In 



(2.24) 



(2.25) 
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Therefore choosing the new scale fj! = Mr with the scale of dimensional transmutation M fixed, the 
renormalization group improvement of the product (2.22) leads to 



yBi\MB\)Vijf{MR) = -5fl(|M«|)Af^ln 

with g = gR{M) , X = Mr/M. 

The expression ( |2.26| ) features a Landau pole at 



Mr, 



M 



g M^x^ In Ixl 
(l-fln|xl) 



(2.26) 



Me 9 



(2.27) 



which, if interpreted in terms of the equation of motion of the scalar field via the effective action (2.20) would 
signal infinite time derivatives when the value of the scalar field reaches the putative Landau pole. Since 
the large N limit does not restrict the coupling to be weak, the value of g can be 0(1). Therefore if the 
dynamical evolution of the expectation value of the scalar field is solely determined by V'g'^j we would expect 
large derivatives and non-analytic behavior of the dynamics as the scalar field approaches the position of the 
Landau pole. 

An important result of this work is that the Landau pole is not relevant for the dynamical evolution of 
the scalar field and contributions from particle production that cannot be captured by the effective potential 
become very important. These non-equilibrium contributions lead to smooth dynamics as the expectation 
value of the scalar field nears the Landau pole. 



3 Large N Yukawa Dynamics 

Having studied the static aspects of the Yukawa theory in the large N limit via the effective potential and 
its renormalization group improvement, we now focus our attention on the dynamical aspects of this model. 
As discussed in the introduction, we consider a system of N fermions coupled to a scalar field 



$. The Lagrangian density in terms of bare fields, mass and couplings is given by equation (2.1) above. In 



order to study the dynamics as an initial value problem we introduce an external "magnetic field" coupled 



to the scalar field so that C C + h{t)^B with C the Lagrangian density in (2.1). The external source h(t) 
serves to generate a spatially homogeneous expectation value for the scalar field 

We will assume that the source was switched on adiabatically in the infinite past and then slowly switched 
off at t = 0. This means that the scalar field evolves in the absence of this external source for t > 0. This 
adiabatic switching on procedure allows us to establish a connection with the effective potential formalism 
of the previous section. If the initial state as t ^ — oo is the vacuum, an adiabatically switched on source 
ensures that the state is the adiabatic vacuum; recall that the (zero temperature) effective potential is 
referred to the expectation value of the field in the vacuum state. We will discuss this issue in greater detail 
below when we address the renormalization aspects. 

In order to extract the dynamics of the mean field, we expand the scalar field as $b {x, t) = ViV 5b (t) + 
X {x, t) with {x {x, t)) = 0. Implementing this last equation within the path integral via the tadpole 
method [p^ , we find that to leading order in 1/N, we arrive at the following equation of motion for the 
mean field 6b (t): 

6b it) + ml 6b it) + -^6% (t) + ^ E = Ht)- (3.1) 

i—l 

Note that this equation is non-perturbatively exact in the N ^ oo limit. The Dirac equations for the N 
species of fermions V'i Eire : 

ii 7*^9^ - Mb it)) = 0, where Mb it) = yB 6b it) ■ (3.2) 
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The above equations are invariant under a permutation of the fermion fields so that we need only deal with 
one of them, denoted generically as V'; therefore we make the replacement 



(V'j ipt) VB ('0V') 



N 



(3.3) 



in eq.(|3.lD. 

In order to proceed, we expand the spinor field operators in terms of a complete set of mode functions 



solutions of the time dependent Dirac equation into eq.(3.2) 



where the creation/annihilation operators obey the usual anticommutation relations 
and the Dirac spinors satisfy the completeness relation 



= 5ab 



with a, 6 being Dirac space indices. Furthermore, eq.(3.2) implies: 

(, - 7 • p - Mb (t)) {t) = 0, {i fdo+i-p- Mb {t)) ^i"' (t) = 0. 



(3.4) 



(3.5) 



(3.6) 



(3.7) 



Since the time derivative operator 9o is singled out by our need to consider the time evolution of the system, 
it is best to work in the basis in which 7" is diagonal. Writing the four-component spinors (t), vi"'^ (t) 



in terms of two-component spinors XiffK^.ip'' j * = I7 2: 



we can use eq.(3.7) to find: 



(z So - Mb it)) h p- (t) Xtj = (<?•?) /2 p- (t) X% 



(a) 



(3.8) 



(3.9) 



(z9o + MB(<))/2p-(t) xg-^ = (^•p)/ip-(t) xfe\ 



(z ao - Mb it)) 51 p- it) eg-' = - (? • 52 p- (i) 
(z ao + Mb it)) 52 p- (<) eg-' = - (^ • 51 p- (0 • 
We also impose the normalization conditions 

cQ/3 _ T.(a)t /.^ T/(/3) 



which together with xg'^ xg' = 5°^^ , i = 1,2, imply that 

\flAtf + \f2pit)f = l 



(3.10) 



(3.11) 



(3.12) 
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which is a consequence of the conservation of probabihty in the Dirac theory, which in turn is a consequence 
of the fact that the Dirac equation is first order in time. This constraint on the mode functions wiU have 
important consequences in the dynamics and underhes the mechanism of Pauh blocking, as will be analyzed 
in detail below. 



We determine the initial state by demanding that ui"'^ {t = 0) represent positive energy states while 



Vp {t = 0) represent negative energy ones: 

t doU^f (t = 0) = Up {t = 0) t/pl") < 0) , (3.13a) 
z 5oFi"' (< = 0) = -LOp {t = 0) y^") (i = 0) , (3.13b) 

where ujp (t) — \J + Mb (i)^ are the mode frequencies. Considering an initial state for which the time 
derivative of the scalar field vanishes at i = (which can always be achieved by a shift in the time variable) 
and then evaluating eqs.(3^,3.10) at < = and using eqs.(3.13), we find: 



(^o-Mo) A,~(0) xif = {S-P) /2p(0) (3.14) 
(c^o + Mo) /2p(0) X^^p^iS- 

where = ^p (0), Afo = Mb (0). This leads to the relations 



(c.0 + Mo) /2p (0) xj^ ^{S-p) /ip-(O) x[% (3.15) 



A similar analysis for the negative energy modes shows that in fact gi p (t) = f^^ (t) , 52 p {t) = flp (t) , ^2^'' = 

((T • p) ^^Z' so that we only need to solve for the positive energy modes. Using the relation between the spinors 
in ( |3.1(]| ), equations ( |3.9[ ) reduce to: 

{ido - Mb (t)) fip{t) = pf2p (3.17) 

iido + MBit))f2p{t)^pflf. 

These can be seperated into two second order equations: 

(92 +p2 + Af2(i) ±,Mb (i)) /i,2 pit) - 0. (3.18) 

The second order equations have twice as many solutions as the first order equations, but since the initial 
conditions for these equations are determined from the first order equations, the correct, physical, pair of 
solutions will be found. The second order equations are of interest because they are more amenable to the 
WKB expansion in the next section. 

The quantity (^'ipip^ appears in the mean field equation eq.( |3.l[ ); we need to calculate it in terms of the 

mode functions {fip{t) , /2p (i)}- This is easily done for our state when we note that ^^^q, &p./3^ — while 
dp,o. 4,a ) = 4/3} ) = (27^)' '^^'^ (P - 9^ use gipit) = /*-(t) , g^pit) = /*-(t): 

d'^p 

Another important quantity for us is the energy density pf in the fermionic fluctuations: 



(V5V') = 2 J P[\f,p(t)\'~\hp{t)\'). (3.19) 



Pf ^ (1^0 ) = \{^H-ijH)^2 j (AV (0 hp it) + r2p it) hp it)) , (3.20) 



(27r) 
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where Tq is the indicated component of the fermionic stress energy tensor operator. 

At this point it is important to highhght the connection with the static effective potential studied in 
section above. 



The equation of motion for the mean field eqn. (3.1 ) suggests the identification of the fermionic contribu- 
tion, the last term of (3.1), with the derivative of the fermionic contribution to the effective potential (divided 
by N) given by eqn. (2.7). Such an identification, however, could only hold for a time independent mean 
field. Indeed in this case, when the scalar field is independent of time, the solutions of the mode equations 
with the initial conditions given in eqn. (3.16) are given by fi.p{t) — /i.p(O); /2,p'(<) — /2.p(0) leading to 



mUiic = -2 / ^f"",,, , (3.21) 

which is the derivative of the fermionic contribution to the effective potential with respect to the effective 
mass Mb- Thus the connection with the effective potential for a static mean field is manifest. As it will be 
explained in detail below, for a dynamically evolving mean field there are non-equilibrium contributions that 
cannot be captured by the static effective potential and that are responsible for a wealth of new phenomena 
associated with the non-equilibrium dynamics. 



4 WKB Expansion and Renormalization 

While the renormalization aspects of the Yukawa theory to leading order in the large N approximation have 
been studied previously ||ll|, here we provide an alternative method, based on the WKB expansion, that 
makes contact with the effective potential formalism studied in section ^ above. 

4.1 WKB Solution of the Mode Equations 

Both the chiral condensate {'tptp) and the energy density pf are divergent. In order to construct the renor- 
malized equations of motion, we need to extract their divergent part s so that they can be absorbed by 



appropriate counterterms. This can be done by finding solutions to eq.( 3.18 ) that allow for a large momen- 
tum expansion, i.e. a WKB expansion. We will concentrate on f2p{t) since the results for fip{t) can be 
then be obtained via the replacement Mb (t) ^ —Mb (t) throughout, including in the initial conditions. 
The WKB ansatz for (t) is: 

(t) ^ A{t) exp i / dun{u) , (4.1) 







where from now on we omit the p index on the mode functions. Insert this into eq.( 3.1§| ) and take the real 
and imaginary parts to find: 

^ + (t) - o2 rn • 2 ^ Mb it) 

-^+u;p[t}-Q (t), 2-^ + ^-^^. (4.2) 
We can solve these equations: 

/f ) (t) = ^=L= exp /* du (tn (u) + ) , where (4.3) 




1 fnit) 

2 [n{t) 

Mb jt) tl (t) 
fl{tf 



+ (4.4) 
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Given f2^^ (i), the second, linearly independent, solution f2^' (t) is 



(in 



it) 

Fit) 



A'^ {t)F{t), 
du 



(4.5) 



i'\u) 



Without loss of generality we can take the lower limit of integration in F (i) to be i = 0, since including 
an arbitrary constant in F [t) will only give rise to a part proportional to ^ (t) in /2^^' (i). This in turn 
can be absorbed when constructing the appropriate linear combination to satisfy the initial conditions. The 
overall factor of —i was chosen so that the Wronskian of f2^^ (t), f2^^ (t) would coincide with its value in 
the equilibrium case. 



Write /2 (t) ~ Aj f2^^ (t) + Ajj /j^^^* (t), and impose the conditions 



to find: 



h it) = /2 (0) 



/2(0) 



A'' it) 

fP (0) 



^0 - Mq 
2^0 



«/2(0) =Cc>o/2(0) 



1+ {Loo + n{o) -i 



2n (0) 



fP (of i^W 



(4.6) 



(4.7) 



4.2 Renormalizing 

We begin the renormalization program by first studying the divergences in the fermion (chiral) condensate 
i^^ip') ■ The ultraviolet divergences of this expectation value can be extracted by performing a high momentum 
expansion of \ f2 it)f — \ fi it)f] the only terms that will be relevant are those proportional tol/fcP,0<p<3 
for large k. An alternative formulation can be found in references l2C|, Furthermore, since we can obtain 
/i (t) from /2 (t) by making the replacement Mb it) ^ — Afg (t), I/2 (i)| ~ |/i (i)| is odd in Mb it). 



We first solve for the WKB frequencies by iterating eq. 
the momentum in the large momentum limit: 

n(°^ (t) = uj (t) 



once and keeping the appropriate power of 



{t)=Uj{t) 1 + 



Mb it) 

4^(^)^ 



(4.8) 
(4.9) 



We will also need the high momentum behavior of F (t) as defined in eq.(4.5). This is most easily obtained 
by integrating by parts a sufficient number of times to extract the relevant part. Doing this for {^p^p), we 
find 




. Ma(0) Mb jO) 
''2n{0)^ 417(0)^ 



(4.10) 



cxp-2 / dt' iVL(t') 



Mb jt') 

2n (f) 



After some algebra we find: 



{\f2it)f-\hit)t 



MBit) , MBit) Mb iO) 



Lu(t) 

Mb (0) 
2ujoUj (t) 



sm 



4^3 (t) AujIuj (t) 



dt'oj {t') 



cos 



0\ 



dt'uo it') 



+ 



(4.11) 
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Note that 



|/2(0)|'-|/i(0)|^ 



|/2(0)|'-|/i(0)|^ 



Mo 



div 



(4.12) 



There are several important aspects of the renormalization of the condensate {4>^) that should be em- 
phasized at this point: 



The momentum integral of the first term in eqn.(4.11) yields the derivative of the effective potential. 
This would be the only contribution in an adiabatic limit in which the derivatives of the expectation 
value of the scalar field vanish. This observation allows us to make a first contact with the preparation 
of the initial value problem via the external source term h{t). Switching this source on adiabatically 
from the infinite past up to the initial time t = leads to the first term in eqn. (4.11) for t < only. 



Integrating the second term in (4.11) in momentum up to an upper momentum cutoff A leads to a con- 
tribution of the form cx 5 sit) ln(A); this should be identified with the wave function renormalization. 



By choosing Mb (0) = we are able to dispense with the fourth term in ( 4.11 ) 



The third term in (4.11), proportional to Mb (0) , has a logarithmic UV divergence at t = which can 
potentially give rise to initial time singularities in the equations of motion pO|, pl| , |2^ , for which there 
are no counter terms. 

However, in the infinite momentum cutoff limit, the contributions that give an initial time singular- 
ity, are actually finite for any t > 0. This is because the integrand is averaged out by the strong 
oscillations At t = the contributions from the second and third terms cancel exactly. For finite 
but large cutoff A it is a straightforward exercise to show that the combination of the second and 
third terms (proportional to M) is finite and small for t < 1/A. The UV logarithmic divergence in the 
combined second and third terms begins to develop on a time scale t > 1/A. Thus, for t ^ 1/A, when 
the third term in (4.11) is finite, we can write 



(^^)^(t) = (^^)>^,(0 + (^^)>(t) 



with ('/' V')(jjv2 (^) given as the momentum integral of the first two terms in cq.( l.ll ) 
Using an upper momentum cutoff A and dropping terms of order (■^) we find 



(^^)div2W = 



27r2 



-M'nit) \n 



Mb. (t) 



2A 



+ - -MRit)A'- 



Mnit) 



In 



Mr [t) 



2A 



(4.13) 



(4.14) 



and (^-tpip) (t) contains the third term of (4.11) and is finite for t ^ 1/A. In section(4.5) we will show 
that this is also finite for t=0 since Mb{0) is proportional to l/ln{A). 
Alternatively we can also write 



{^^)B = {^^),i.At) + {^'^)AF(t) (4-15) 



with (^ipip'^^^^^ (t) given as the momentum integral of only the first term in eq.(4.11) 



^This can be seen simply by considering the contribution to the momentum integral in the limit of fe >> Mg. The 
contribution of the large fc-modes can be estimated by taking the ultraviolet cutoff to infinity but introducing a lower momentum 
cutoff ^, leading to cos(2kt)^ = — Ci(2/it) with Ci{x) the cosine integral function which is finite for t > and diverges 
logarithmically as t — * 0. 
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Mr (t) 



2A 



(4.16) 



Now (^ip'4')j^p (t) includes both the second and third term in eq.( l.ll| ). It is clear from the discussion 
above, that (i/i^)^^ (t) is actually finite for t < 1/A; it vanishes identically at t = and does not 
have any initial time singularity. This quantity will, however, develop a logarithmic divergence due to 



the second term in eq.(4.11) associated with wave function renormalization for t ^ 1/A (hence AF for 
"almost finite"). 

Several important features of the above expressions must be highlighted. First, the argument of the 
logarithms contains the full time dependent mass, unlike a renormalization scheme that only extracts 
the logarithmic divergences in terms of the initial mass([0 |l^, |2l|, ^); note that these schemes differ 
only by finite terms. However, as will become clear below, keeping the full time dependent frequencies 
in the denominators will lead to the instantaneous effective potential, i.e, the static effective potential 
but now as a function of the time dependent mean field. Furthermore, taking the full time dependent 
frequencies will lead to the identification of the effective coupling that "runs" with the amplitude of the 
mean field. This identification will allow us to establish a direct correspondence with the RG improved 
effective potential. In this manner we will be able to clearly separate the contribution from an adiabatic 
or instantaneous generalization of the static effective potential, from important non-equilibrium and 
fully dynamical effects which can only be described in terms of the time dependent mean field. In 
particular, we seek to clearly separate the effect of particle production and its concomitant contribution 
to the dynamical evolution. 



Second, the term tt-j- 



-Ml 



{t){ 



In 



2A 



+ i) - Mnit)A^\ in both {i^^)^^^ (t) and (^V')^;,! W 
will lead to the effective potential. This is implied because this term does not depend on the derivatives 

of MBit). 

Finally, while {\}}ipYfi{t) is finite for t 3> 1/A but features an initial time singularity at t = 0, 
(xp tp) [t] vanishes identically a,t t — Q, it is finite for t < 1/A but exhibits a logarithmic ultra- 
violet divergence cx (51n(A) associated with wave function renormalization for t 3> 1/A. 

If we insist in using the split ( |4.13| ) , and thus extract the wave function renormalization divergent term 
at all times, including t < 1/A the quantity (iptp) ^ wi ll contain an initial time singularity given by the 
short time limit t < 1/A of the third term in eq.(4.11). 

In references [^0|, ^ ^ this initial time divergence is dealt with in several ways: by choosing the 
initial conditions on the mode functions to a high (fourth) order in the adiabatic expansion or by 
performing an appropriately chosen Bogoliubov transformation of the initial state pl|, or equivalently, 
by including a counterterm in the external "magnetic field" h (t) so as to cancel this singularity at 
t = 0. All of these methods are equivalent and lead to a set of equations that conserve energy and are 
free of initial time singularities |pO|, pT| , 

However, these methods all suffer from the drawback that they do not lead to an interpretation of the 
dynamics in terms of the effective potential. This is evident in the Bogoliubov approach advocated 
in references [^l[ ^ since the Bogoliubov transformation involves the derivatives of the mean field, 
and the Bogoliubov coefficients multiply the terms that lead to the effective potential, thereby mixing 
terms that depend on the derivatives of the mean field with terms that arise from the adiabatic effective 
potential. We refer the reader to reference]^ for a thorough exposition of the Bogoliubov method. As 
discussed in detail in|^, |2^, ^ any approach to regulating the initial time singularity of (^'pip'^^ (t) 
(when extrapolated to i = 0) is tantamount to a redefinition of the state at i = 0. 

Instead of seeking a regularization of this initial time singularity, we recognize that this singularity 
arises from trying to extract a wave function renormalization from very early time even when there is 
no such divergence. We interpret the fact that the logarithmic divergence associated with wave function 
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renormalization emerges at time scales i 3> 1/A as the build up of the wave function renormahzation 
over this time scale. This is consistent with the adiabatic hypothesis of preparation of the state at 
t = via an external current. At this point the following question could be raised: in the formulation 
presented above, where is the adiabatic assumption explicit? The answer to this important question 



is in the initial conditions on the fermionic mode functions given by eqn.(4.6)(up to an overall trivial 
phase). These initial conditions determine that the state at f = is the vacuum state for free Dirac 
spinors of mass Mb{0) which is determined by the value of the mean field at i = 0. Obviously this 
is the state obtained by adiabatically displacing the mean field from the trivial vacuum. As we will 
see in detail below, recognizing that {{pi') j^p {t) vanishes identically at i = and is finite over a time 
interval i < 1/A will allow us to solve the initial value problem self-consistently with the result that 
the potential initial time singularity is simply not there and the evolution is continuous throughout. 

Thus in summary, the above discussion highlights the very important dynamical aspect that the wave 



function renormalization builds up on time scales t > 1/A and eqn.(4.13) or (4.15) must be used according to 
the time scale studied in the evolution. Ho wever , as will be discussed in detail below, when we study the full 
equations of motion, we will see that eqn. ( 4.15| ) is far more convenient for numerical studies. Furthermore 



we will find in section(4.5) below that a self-consistent analysis reveals that in fact there is no initial time 



singularity, i.e, (ipip) ^ (0) is actually finite 



4.3 The Renormalized Equations of Motion 

We now have all the ingredients to obtain the renormalized equations of motion. We recall our system of 
equations: 



^ + p2 + Mlit) ± I Mb (t) ] h.2 P (t) 
Sb (t) + ml 5b (t) + If 4 it) + VB = 



(4.17) 
(4.18) 



where we have set h{t > 0) = 0. 

Starting with the fermionic mode equations, we see that since there are no operators present to absorb 
the divergences in Mb (t), we need to impose the condition: Mb (t) = Mr (t) or ysfe (t) = yRSn (t). When 
coupled with the wavefunction renormalization conditions below, this will relate the bare and renormalized 
Yukawa couplings. This condition also implies LUp (t) ^ — tOp {t) p . 

To study the dynamics for time scales i 3> 1/A we use eqn. ( [4. 13 ) to separate the divergences, leading to 
the following equation for the mean field 



4^2 




-0, 



(4.19) 



where /i is a renormalization point. Set Sb (t) — VZSr (t), so that AIb (t) — Mr (t) implies yB — yR/VZ, 
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and choose the coefficient of Sr to be unity. This yields the renormahzation conditions: 

2A^ 



1+ Ai^^ 



= 1, 



Z 



27r2 



Ab 
3! 



4ln^ 
27r2 fi 



3! 



Zys = VR (m) 



(4.20a) 

(4.20b) 

(4.20c) 
(4.20d) 



These are exactly the same renormahzation conditions obtained from the renormahzation of the static 
effective potential, together with the wave function renormahzation obtained from the one loop fermionic 
self energy calculated in section (^ above. Furthermore, the renormahzation of the Yukawa coupling (4.20d) 
guarantees that the fermionic mode equations are renormalization group invariant since they depend only 
on the product ysSs = yR^R- 

The renormalized equation of motion for the mean field now becomes 



1 - 



9r M 



9R it^) 
2 

Afl (a*) 



127r2 5|(^) 



In 



In 



cMr {t) 



Mr 



.{t^)MR^ 



(4.21) 



Mr (t) 


4) 


A* 





A4(t) + 2^2^^^^> 



0, 



where gR(ij) = y'jii^) /2vr2, as defined previously, and we have multiplied through by yR to write the 
equation of motion for the renormalization group invariant product yRSR. We can rewrite this as: 



1 



/2m?, 



(m) 



D{ti) V 5i? (Ai) 



1 Afl(A^) 

12^2 g2^ (^) 



- In 



MR{t) 


-I) 
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,(<) + 47r2(v;,A)^ =0, (4.22) 



where D (fi) = {2/gR (^) - In {cMr (t) 
eq.(4.22) are ji independent: 



Using the following relations we can see that all the terms in 



TO 



9R (m) 
2 



B 
9B 



9R (m) 
Afl (m) 



2 , 2A 
= — + In — 
9B M 

1 Af 



12^2 



In particular, 



9R (m) 



In 



12^2 ^1 



cMr (t) 



In- 



2A 



<?fl(eMfl(t)) 



(4.23a) 
(4.23b) 
(4.23c) 

(4.24) 



Since gR (/Lt) oc yj^ (/i) > 0, unitarity appears to require that the range of validity of the theory be restricted 
to be below the Landau pole at El = e~^/iexp {2/gR (/i)). As discussed in the Introduction, however, this 
may not actually be necessary under all circumstances; we will explore this issue further in the next section. 

We can now compa re to the static case. Taking ttir = (as in the static case) we see that the term 
proportional to in ( 4.22 ) is precisely the derivative of the static effective potential given by eqn. ( 2.15 ) 
but now in terms of the time dependent renormalized mean field. Therefore, just as with the static effective 
potential, we introduce the dimensional transmutation scale M by demanding that the instantaneous or 
adiabatic effective potential, i.e, the static effective potential in terms of the time dependent mean field, have 
a maximum at this scale when tur — 0: 



1 Afl(A^) 
127r2 5|(^) 



, M 1 

-In = 0. 

[i 2 



(4.25) 
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We also define x = Af^ (t) /M, g = gj^^M) ,q = k/M, ni^ ^■m\{ijL^ M) /NP, t = Mt. Then eq.( |4.22| ) can 
be written as: 



X" (r) + 9R (ex (r)) (r) - (r) In |x (r) | + 2^^ ^^^|f^ j = 0, (4.26) 

where primes denote r derivatives. 

For mn = we see immediately that the combination — .'7j?.(ex(' '"))x'^ (''") In lx (''')l is the derivative of the 
renormalization group improved effective potential (see equation ( 2.26| )) in terms of the running coupling 
constant at a scale x{t) (np to a finite term). 

Thus this form of the mean field equation of motion splits off the effects due to the RG improved effective 
potential ( the first two terms) from those due to the time evolution of the fermionic fluctuations which 
are represented by (^V-')j^- We also see that the effects of wavefunction renormalization are encoded in the 
prefactor gji (ex (t)) multiplying the potential and fluctuation terms. 



A remarkable aspect of equation (4.26) is that the effective coupling depends on time. In a well defined 
sense, this is a dynamical renormalization much in the same way as that explored in real time in references [p3[. 
As mentioned before we could have renormalized simply by absorbing the divergences with the frequencies 
at the initial time. However, doing this would still leave "large logarithms" arising from when the amplitude 
of the mean field becomes large. As the mean field evolves in time, the fermion fields probe different energy 
scales. Since the coupling "runs" with energy scale, it is then natural that it will "run" as a parametric 



function of time through the evolution of the mean field. This physical picture is manifest in eq.( 4.26 ). This 
important aspect of our study is a novel result which only becomes manifest in a real time non-equilibrium 
framework that allows one to study the dynamics of the fully renormalized fields and couplings. 

It would appear from this equation that the dynamics of x could be singular asx crosses E^/M since 

where 



E=^=e~^~'^. (4.28) 
M ^ ' 

This would certainly be the case if the equation of motion only involved the derivative of the static effective 
potential. However, we will see below that this is not the case when the fluctuations are taken into account. 

We conclude this section by collecting together our renormalized equations of motion and initial con- 
ditions, now in terms of dimensionless variables, valid in principle for t ^ 1/A when {ipip)^{T) is free of 
potential initial time singularities: 

^ + q^+X^r)±i x' (r)) /i,2 ? (r) = 0, (4.29) 

x"{T)+9R{ex{r)) (^^X (r) - x' (r) In |x (r)| + (r)^ = 0, (4.30) 

_o_2(MiM_o r^„„2 A, .m2 ,.,2 / X(r) , X"ir) 



with the initial conditions 



/i,2,-(0)^y '^^^"]^'^Q^^"^ , */i,2?(0) = c.,(0)/i,2,-(0) (4.32) 

X(0) = X0 , X(0) = (4.33) 
(^V')^^(0)-0. (4.34) 
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The quantity T{t) is the fermionic condensate ('0V') after subtracting the contributions that lead to 
the effective potential and the wave function renormalization. It therefore represents the purely dynamical 
fluctuations in the fermion field associated with particle production. 



4.4 RG Improved Effective Potential 

Since we will focus our study on the effects of the dynamics, we start by highlighting the dynamics that 
would ensue from a consideration of the RG improved effective potential aZone|2|. A comparison between 
the renormalization group improved effective potential eq. ( 2.26| ) and the renormalized equation of motion 
(4.30) clearly indicates that we can extract the dynamics that would ensue solely from the effective potential 
by neglecting the dynamical contribution of the fermionic fluctuations encoded in J-{t) in the equation of 
motion (4.3C) leading to 



X"(r) 



X^(r)ln|x(r)| =0 



(4.35) 



2/g-ln|ex (r)| 

Notice that the denominator in the second term comes solely from wavefunction renormalization consider- 
ations and that it has a zero at the Landau pole, as discussed above. This truncated form of the equation 
of motion displays clearly the connection with the renormalization group improved effective potential, as 
obtained in section (2.1), eqn. (2.26)). Obviously for generic values of fh? and the numerator will not 
vanish when the denominator does, leading to infinite acceleration at the time when the mean field reaches 
the Landau pole. Thus, an investigation of the above equation would lead to unphysical behavior of the 
mean field as it approaches the Landau pole and one would then conclude that the existence of a Landau 
pole precludes a sensible interpretation of the theory when the amplitude of the mean field is comparable to 
the position of the Landau pole. 

We bring this discussion to the fore because it is one of the important points of this study that the effect 
of the fluctuations is very dramatic and completely changes the picture extracted from the effective potential. 



4.5 Full dynamics: fermionic fluctuations 

Having established the connection with the re norm alization group improved effective potential, we now study 
the evolution of the full equation of motion ( 4.30 ) including the fermionic fluctuations T {t). 

However at this point we face two problems: i) the renormalized equation of motion ( 4.30| ) is only valid 
for t ^ 1/A and cannot be extrapolated to the initial time t = Q because of the potential initial time 
singularity in discussed above. This in turn, entails a potential problem with the initialization of 

the dynamical evolution, ii) Equation (4.30) cannot be numerically implemented accurately because !F{t) 
requires a subtraction that involves the second derivative of the mean field at the same time as the update 
(see eqn. ( 4.31 ). Both problems can be circumvented at once by invoking the split ( 4.15 ) or equivalently by 
introducing the dimensionless quantity J^af{t) as 



T{t)=TAv (r) 



^AF (r) - 2 




(4.36) 
(4.37) 



We call Taf (t) the "almost finite" fiuctuation, since it only has the part of the divergences proportional 
to X {t) subtracted. Furthermore, as discussed above, J^af (0) = 0. 

We emphasize that we have not changed the equation but have merely redistributed the terms, just as 
the two alternative forms of writing the fermionic condensate given by dil^ ) and ( |4l5| ). 

The integral in eq.( 4.36| ) can be done explicitly. We can then combine the terms proportional to x" {t) 
and rewrite eq. ( 4.3C| ) as 
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X"(r) + 



Za (t) = 2/g-l + 



■x(T)-x'(T)ln|x(T)|+^AF(r) ) -0 
A 



(4.38) 



V^'+X' (r) 



In (a + v/A2 + (r)) 



We have kept the full expressions for the integral for the sake of completeness. However, if x (t) approaches 
the cutoff, fermion modes with momenta of the order of the cutoff will be excited indicating that we are 
approaching the limit of validity of our numerical approximation. This means that we should restrict x (t) 
to be much less than A. Doing so allows us to simplify the expression for Za (r): 



Za{t) -2/5-ln2A = ln 



E 



2e-iA 



In- 



E 
A' 



(4.39) 



The important point to note here is that Za (t) is independent of time. In particular, in this formulation 
the Landau pole does not give rise to any singular behavior. 

From the renormalization of the coupling g given by eqn. ( 4.23b| ) we identify 



2 



9B 



(4.40) 



i.e, the coupling at the scale A or the "bare coupling" . 

We reiterate that we are still solving the renormalized equations; all we have done is reformulate them 
for the purposes of numerical analysis. The fact that Za (r) is non-vanishing provides a strong hint that 
there should be no problems as x crosses the Landau pole. We shall see below that indeed this is borne out 
by the numerics. 



The equation of motion (4.38) is now in a form that can be studied numerically. In particular, it can be 
initialized for any large but finite cutoff. Furthermore, the form of the equation of motion given by eqn. ( 4.38| ) 
above suggests that the dynamics is smooth even at time scales f S> 1 /A when J-af {t) develops an ultraviolet 
logarithmic divergence. This is so because the logarithmic divergence in J-'Api''') will be compensated by 
the logarithm in Za in the denominator. Hence we conclude that the equation of motion ( 4.381 ) has a well 
defined initialization and the dynamics is smooth without logarithmic divergences or discontinuities as x(''") 
approaches and passes the Landau pole. This in fact will become clear from the detailed numerical analysis 
provided in the next section. The equivalence of the equations of motion in terms of J^af or and the 



observation that the equation of motion (4.38) does not feature any initial divergence suggests that the 
equation of motion (4.30) is free from initial time singularities and the numerical evolution is indeed smooth. 
That this is indeed the case can be seen as follows. From the fact that !Faf (0) = we now find that 



2 / rh^ 

X (0) = -^ — x(0)-x' (0)ln|x(0)| 
Za \ 9 



Furthermore, from the relation between and J-^af given by eqn. (4.36) we find that 



m^^xio) 



1 - In 



2A 



In 



In 



9 



-X(0)-X' (0)ln|x(0)| 



(4.41) 



.42) 



In the large cutoff limit, the above expression becomes cutoff independent and we conclude that the equation 
of motion in terms of the renormalized fermionic condensate is also free of initial time divergences. Thus 
either formulation can now be used for a numerical study with well defined initialization and smooth evolution 
throughout with no cutoff dependence in the limit when the cutoff is taken to infinity. 

The resolution of the initial time singularity is now clear: from equation (4.11), it is clear that the initial 
time singularity is completely determined by x"(0) which in turn must be obtained self- consistently from the 
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renormalized equation of motion in terms of T . The formulation of the equation of motion in terms of T af 
which vanishes at t = allows us to find the value of the second derivative at the origin, the logarithmic 
singularity is now encoded in Za which leads to an initial value of the second derivative of the mean field 
that is vanishingly small in the limit of large cutoff. 



The solution (4.42) has a remarkable aspect that explains how the dynamics manages to be smooth when 



the amplitude of the mean field approaches the Landau pole. Consider the initial value problem in which 



x(0) is very near the Landau pole, i.e, x(0) ^ this case we find from (4.42) that 



no) = - ( — X (0) - (0) In Ix (0) I ) (4.43) 



which leads to the conclusion that the coefficient of the coupling 5fl(ex(0)) in ( 4.3C ) actually vanishes. That 
is to say, the potential singularity at the Landau pole is actually rendered finite by an exact cancellation 
between the fermionic fluctuation and the adiabatic effective potential. We will see numerically below that 
this remarkable feature is borne out by the dynamics in all cases. 



In the cases below we can solve the mean field equation in the form given by eq.(4.38) together with 
the fermionic mode equations, simply because the update does not require the specification of the second 
derivative and is therefore more accurate. However after each step in the iteration we have constructed !F 
and checked that the value of the second derivative obtained from both formulations coincide, thus providing 
a numerical check of consistency. 

Standard numerical techniques (fourth order predictor-corrector Runge-Kutta ODE solver together with 
a fourth order Simpson's rule integrator) are used. We also compute the fermion occupation number relative 
to the initial vacuum state in each momentum mode as a function of time given by[T]|| 

n, (r) - I /2 (r) f*^^ (0) - /i (r) (0) | ' (4.44) 

In the section (^) below, we analyze the equations governing the system for various values of the param- 
eters rh'^,g,xo- 

4.6 Renormalizing gf 

Before proceeding to the numerical study of the renormaliz ed eq uations of motion we now turn our attention 
to the renormalization of the energy density gf. From eq.( 3.20| ), we see that a time derivative of the mode 



functions is involved in computing Qf. This has the effect of bringing down one more power of momentum 
into the integrand and implies that we need the WKB expansion of the mode functions to order k~'^. This 



entails both one more iteration of the WKB frequency, eq.(4.4), as well as one more integration by parts on 
the function F {t). Doing this yields 



(AVW Ap(t) + /2V(*)/2pW), =-u;it)-^^^^^-^^+ (4.45) 

^ ^ div {^t) 



Mb jt) 
Aw [tf 



Mb (t) - Mb (0) cos 2 dt'uj {t')^ - 

Mb (0) (Mb (<) - Mb (0)) rt 

^ ^■sin2 / dt'uj{t'). 



8lo {tf Jo 



When the factor of two in eq.( |3.2C ) is included, the first term above gives the contribution of the zero 



point energy densities for a four component Dirac fermion. This would be the one-loop approximation to 
the effective potential Q if the frequencies were constant. In this case, we can consider this piece to be an 
adiabatic approximation to the effective potential as described above. The other terms have been grouped so 
that their contribution vanishes at the initial time. While the fully renormalized energy density has a lengthy 
but not very illuminating expression, we just highlight its main features. After mass, coupling, wavefunction 
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renormalization and a substraction of the zero point energy (time independent and proportional to the fourth 
power of the cutoff) the energy density is finite and conserved by use of the fully renormalized equations 
of motion for the scalar field and the mode functions. There is no unambiguous separation between the 
fermionic and scalar energy density because of the wave function renormalization, which arises from the 
fermionic fluctuations but contributes to the kinetic energy of the scalar field. The total energy density 
is renormalization group invariant, finite and conserved. Furthermore we have checked numerically in all 
cases that the energy density is constant throughout the evolution to the accuracy required in the numerical 
implementation, thus providing an alternative check of the reliability of the numerical calculation. 



5 Solving the Equations of Motion 



We now turn to a discussion of the actual time evolution of the coupled scalar-fermion system. At this stage 
we summarize the discussion of the previous section on renormalization to be able to focus on the important 
aspects to be gleaned from the numerical study: 



The fully renormalized equa tions of mot ion fo r the fermion field modes ( 4.29| ) and for the mean field 
either in the form given by (4.30) or by (4.38) are free of initial time singularities or ultraviolet diver- 
gences. They can be consistently initialized and lead to smooth evolution. The two forms of the mean 
field equation are completely equivalent as they are obtained one from the other by a rearrangement 
of terms. While eqn. (4.3C) seems to suggest a singular behavior when x reaches the Landau pole, the 
equivalent form ( 4.3^ ) suggests smooth and continuous evolution. 



• Mass, coupling and wave function renormalization and a renormalization of the zero point energy 
(time independent) renders the energy density finite and conserved as a consequence of the equations 
of motion. 



We now study in detail different cases to bring the role of the fermionic fluctuations to the fore. 
5.1 Case 1: ^ , A < E and |xo| < Xmax 

For this case we will use rh? / g — l,g — 0.261 and look at two different values of xo • Xo — 0.75,0.2. For 
reference, the location of the maxima of the effective potential are at |xmax| = 1-53 so for these values of xo 
we expect x (t) to oscillate about the origin. Figure (^) below displays xi^) ^ ^ function of r for xo = 0.75. 

The first figure reveals that the amplitude of the mean field decays as would be expected; there is energy 
transfer to fermion particle production. These two figures taken together exhibit an interesting feature. There 
is a remarkable similarity at later times between the mean field and the oscillations in J^. A comparison of 
the two figures suggests that the amplitude of the fluctuations is proportional to %. Since the mean field 
has a decreasing envelope while the fluctuations increase, if such a proportionality exists, it must involve a 
coefficient that is slowly increasing in time. We can actually extract more information from a parametric 
plot of x" versus x, which is shown in figure (^. 

The fact that such a tight curve is produced is indicative of an underlying relationship. In fact, we found 
that this curve could be well fit to ax (t) + bx^ (t) In |x (t) | with a; b very slowly varying over the time scale 
of the oscillations. This in turn implies that the fermionic fluctuations J- (t) can also be fit to this form, with 
coefRcients that are slowly varying function of time indicating that the growth of the fermionic fluctuations 
results in a time dependent correction to the mass term and quartic coupling of the mean field. This, we 
believe, is a noteworthy aspect of the dynamics: the non-equilibrium fermionic fluctuations, those that were 
not accounted for by the adiabatic effective potential, introduce a slow time dependent renormalization of 
the parameters of the effective potential, mass and quartic coupling. 

There is a new time scale emerging from the dynamics that is associated with the (slow) time evolution 
of this renormalization and the decay of the mean field. A full analysis of these time scales is beyond the 
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Figure 3: x{r) and {^^P{T))> vs. r. x(0) = 0.75 
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Figure 4: %"(''') vs. x(''') 



scope of this article, but we expect to use the methods of dynamical renormalization group p3[ to investigate 
the relaxation of the mean field in future work. 

We next consider the behavior of the fermion occupation numbers as a function of time. Figure (|) below 
gives three snapshots of n (k) versus k at different times. 

What should be apparent from these pictures is the existence of band structure in the fermion occupation 
numbers. In fact, the produced fermions seem to have wavenumbers that lie within a region spanning A: = to 
k = kc where k^ xo- While we are used to band structure in bosonic theories with parametric resonance, it 
is unexpected to encounter this structure in fermionic theories. Pauli blocking prevents exponential particle 
production since the occupation numbers can at most become unity. In section (|^) we study in detail the 
issue of parametric amplification versus Pauli blocking. 
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Figure 5: ri(fc) vs. k at times r = 15.01; 51.05; 99.09 respectively. 



Now consider the case where xo = 0.20. 

Figure (^) shows the phenomenon of "catalyzed regeneration" or "revival" first observed in [Q. The 
mean field decays for a while and then regenerates itself. It is important to note that it never regenerates back 
to the original value, always to something less than that. The reason for this can be seen from the evolution 
of the momentum distribution below (Figure (^) which clearly shows an almost saturated distribution of 
particles for momenta A: < xo- 

The allowed band fills up to saturation very early on. After that, the energy in the scalar cannot be 
transferred efficiently to fermions and in fact the fermions, which only couple to the mean field in our 
approximation, begin to transfer their energy back to x- This depletes the band but not completely, which 
accounts for incompleteness of the regeneration. 

5.2 Case 2: m = , A < i? and |xo| < Xmax 

The massless case leads to a qualitatively similar dynamics of the mean field and the fcrmionic fluctuations 
but with definite quantitative differences in the time scale of oscillations damping of the mean field and 
growth of the fcrmion fluctuation as compared to the massive case. These are shown in fig.(^ below. 

We note that now the only scale in the problem is completely determined by the initial value of the mean 
field x(0) whereas in the previous case there were two scales. 

The figures for the dynamics of the mean field and fermion fluctuations are qualitatively similar to those 
of the previous case. Figure (j^) shows the momentum space distribution of fermions. 

The band structure is still present and is still set by the initial value xo- This is as it should be , especially 
since xo is now the only scale in the problem. Finally we again plotted x" versus x with a result identical as 
that shown in figure We can again find a good fit to ax (t) + hy^ (t) In |x (t)| with the parameters a, h 
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Figure 6: x{t) vs. t for x(0) = 0.2. 
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Figure 7: n{k) vs. k at times r = 15.01; 51.05; 99.09 respectively. 
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Figure 8: x{r) and {■l|J^P{T))> vs. r. x(0) = 0.75 




very slowly varying on the time scale of the oscillations. This fit shows that a time dependent mass has been 
generated by the dynamics. This is not so surprising as scalar masses are not protected against radiative 
corrections. 

Again, just as in the previous case we find that new time scales associated with damping in the amplitude 
of the mean field and the renormalization of the parameters emerges. 

5.3 Case 3: fh = , A < E and |xo| > Xmax 

In this case the mean field begins to the right of the maximum of the effective potential and rolls down, 
running away up to the scale of the cutoff, with a similar behavior for {^^{t)). This behavior is displayed 
in figure (^0|). Once the mean field reaches an amplitude of the order of the cutoff scale there is particle 
production on this scale so the momentum integral over the fermionic fluctuations is no longer accurate. 
Therefore the dynamics is no longer trustworthy and the evolution must be stopped. In this case the 
dynamics is what would be expected from the effective potential description: the mean field runs away as a 
consequence of the effective potential being unbounded from below. 
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Figure 10: x(t) and {'4)4'{t))^ vs. t for x(0) = 1.1. 



5.4 Case 4: 



, A > i? and 



In the previous section we have discussed the possibility that the presence of a Landau pole in the fully 
renormalized equations of motion might not be a signal of discontinuities or singularities in the evolution 
and that the dynamics will be smooth when the mean field reaches and crosses the Landau pole. There are 
two important hints that led us to this con clusio n: a:) the Landau pole does not appear explicitly in the 
equation of motion written in terms of !Faf (L38) (the logarithmic divergence in !Fa f is c ompensated by the 
logarithm in Z\). b:) The self-consistent solution for x"(0) and for !F{0) given by ( 4.43 ) when x(0) is very 
near the Landau pole, shows manifestly that the fermionic fluctuations cancel exactly the contribution from 
the effective potential thus making the "residue" at the Landau pole vanish exactly. While this remarkable 
cancellation has been gleaned in a particular case, that in which the initial condition places the amplitude of 
the mean field at the Landau pole, the combination of both arguments are suggestive enough to conjecture 
that the dynamics will be smooth in all cases as xi''') crosses the Landau pole. 

In order to probe this conjecture, we need to consider the case A ^ _E so that x can can evolve past the 
Landau pole but the dynamics should be reliable in such a way that the amplitude of the mean field must 
always be much smaller than the cutoff. 

For A ^ i? we see from the formulation (4.38) that Za < thus for very early time when !Faf ~ 
the acceleration is negative and the mean field climbs the potential hill instead of rolling down. This a 
consequence of the fact that for a theory with a Landau pole, if the cutoff is taken much larger than the 
Landau pole, the "bare" coupling becomes negative but very small. Thus in the "almost finite" formulation 
the only hint of the presence of a Landau pole is in the opposite sign of the second derivative of the mean 
field as compared to the case in which the cutoff is well below the Landau pole. 

Thus initially we expect that the second derivative is very small and negative, the mean field will slowly 
climb up the potential hill and the fermionic fluctuations will grow. 

In what follows we study the cases with g = 1.0; xo = 1.2, A = 200; El ~ 2.7 and g = 0.5; xo = 1.2, A = 
500; El = 20.1 respectively, to illustrate the main features of the dynamics. 

The full dynamics for the mean field and the renormalized chiral condensate in these cases is displayed 
in figures ( [ll| , |l2|) for the mean field and the fully renormalized fermion fluctuation {'4''4>{t))'^- 

These figures show a remarkable behavior: x initially climbs up the potential, reaches the maximum and 
begins falling down towards the mininum, overshoots, climbing up to the maximum on the negative side and 
finally begins a plunge down the potential hill on the negative side. As the amplitude of the mean field reaches 
the Landau pole at a time tlp the (fully renormalized) fermion condensate J^{t) exactly compensates the 
contribution of the effective potential thus cancelling the singularity of the running coupling at the position 
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Figure 11: x(t) and {tjjil){T))> vs. r for x(0) = 1.2; g = 1.0; = 2.7 
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Figure 12: x(t) and (■(/''0(r))> vs. r for x(0) = 1.2; g = 0.5; £'l = 20.1 



of the Landau pole. 

When the amphtude becomes larger than El the mean field begins an oscillatory motion about the 
origin with greater and greater amplitude. The key to understanding this behavior lies in the figure ( p^ ) 
which depicts x" vs x parametrically. Note that the two different cases we are examining are almost 
indistinguishable. 

A numerical fit to an excellent accuracy reveals that for r >> r^p 



X"(r) ^ hx^T) 



(5.1) 



with b slowly varying in time and saturating at long times at a value h < —2.0. 

Thus for T » tlp the fermion condensate provides a small renormalization of the coefficient of x^ In |x|- 
Therefore when the amplitude of the mean field is much larger than El, the J- (r) can be neglected, the logs 
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Figure 13: x" vs. x parametrically for both cases. 



cancel, and the equation of motion takes the form: 

X" ir) + ^T^nf^ (-X^ (-) In IX ir) \+Hr))^0 or 

x"(T)~fex'(T) ; with 6 < -2.0 (5.2) 

This shows that the mean field behaves at late time as if it were in a quartic potential; this would also 
be true in the massive case since the mass term would eventually become subdominant. The "potential" 
appears to "open up" with time i.e. the coefficient of the quartic term becomes smaller and smaller as the 
log in the denominator dominates until it saturates. 

Suppose that the initial value of the mean field were near the maximum and that the coupling is such 



that El >> 1- In this case from the self-consistent solution (4.42) we see that T ~ C(l) (for A >> El). 
However we have argued and seen explicitly numerically that as the mean field approaches and crosses the 
Landau pole, the fermion condensate J^(t) cancels the term x'^(''') In |x(''')l-This implies that the fermion 
condensate becomes of 0{E\). Since the renormalized fermionic condensate is obtained by subtracting the 
terms from the effective potential and the wave function rcnormalization, the only manner in which this 
can actually happen, roughly speaking, is that the mode functions /i,2 evolve in time to almost saturation, 
i.e, that |/i,2p ~ 1 for wave vectors up to the Landau pole scale '-^ e^ls~^ (or beyond). This will lead to 
a fermion condensate of order E\. This argument is shown to be correct by figure (|lj) which displays the 
occupation numbers of the produced fermions as a function of momentum for large r when the mean field 
oscillates. 

Comparing the width of the band of occupied wavevectors with the amplitude of the mean field in the 
oscillatory phase (see figures ^ 12) reveals that the width of the band is proportional to the amplitude of 



the mean field, which, in turn is larger than the position of the Landau pole. Thus states with wavevectors 
up to El (or somewhat larger as shown by fig.(|l^)) are almost saturated with occupation one (per spin). 
We interpret this as the formation of a very dense fermionic plasma with a "chemical potential" of the order 
of the Landau pole, since all states up to this scale are filled by produced particles. 

Summary of the numerical analysis: we want to highlight several noteworthy features that emerge 
from the numerical analysis: a: as presaged by the discussion above, the mean field evolution is completely 
non-singular, even as it crosses the Landau pole. The avoidance of singular behavior is entirely due to 
the dynamical growth of the fermionic fluctuations: they ensure that the quantity —x" (t) In \x {t) \ + T (r) 
vanishes when 2/5 — In | ex (r)| does, b: The backreaction of the fermions provide a small rcnormalization 
to the effective potential after the mean field crosses the Landau pole, and as a result, the true effective 
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Figure 14: n{k) vs k. Left figure is for t — 1.89, g = 1.0, El = 2.7, the right figure for r — 0.35, g — 0.5, i?L = 
20.1. 



potential at large values of x is upright and quartic. The equation of motion is then of the approximate 
form x" ~ ~2x^ where the factor —2 receives small corrections from the fermion backreaction. c: The 
non-equilibrium evolution results in profuse particle production. The fermion occupation number saturates 
generally up to momenta larger than the Landau pole, which in turn results in a very dense medium akin to 
a fermi gas with chemical potential of the order of the Landau pole. 

6 Pauli blocking vs. parametric amplification: 

An important question that we address in this section is the role of parametric amplification in a fermionic 
theory. 

In a bosonic theory, when the expectation value of the scalar field oscillates around a minimum the 
fluctuation mode functions obey an equation of motion with an effective oscillatory mass squared term. 
The resulting Schrodinger-like equation for the mode functions results in a spectrum with forbidden bands 
with positions that depend on wave vector k of the mode function and the specific details of the oscillatory 
potential. For generic initial conditions, the wavevectors k in the forbidden bands lead to an exponential 
growth of the mode functions with an exponent, given by the imaginary part of the Floquet index, that 
depends on the details of the potential. This is the phenomenon of parametric amplification. The exponential 
growth of the mode functions is associated with the build up of a non-equilibrium, time dependent distribution 
function for the bosons, i.e, particle production. 

In the fermionic case, the Pauli exclusion principle restricts the maximum occupation number of fermion 
states for a given wave vector to be one. This is Pauli blocking, i.e, the quantum states must have a 
finite occupation number. Obviously, even if the fermionic fluctuation mode equations contain an oscillatory 
time dependent mass, Pauli blocking should prevent any parametric amplification and forbidden bands from 
occurring. This is somewhat problematic, however. As our detailed numerical work shows, there are cases 
in which the mean field is oscillatory, i.e, x(t) ^ xo cos(fit) and the fermion mode equations ( ^.18 ) take the 
form 



(f 

-^+P^ +xl cos^ {nt) T i^Xo sm{nt) 



fiMt) = (6.1) 



At first glance it would appear that the oscillatory terms would drive parametric amplification of the 
mode functions, which would lead to an unbounded distribution function, contrary to Pauli blocking. This 
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is not the case, though. We now provide an rigorous proof that ferniionic mode functions do not have 
forbidden band structure hence there is no parametric amphfication, and a perturbative argument that not 
only confirms the exact proof but also illuminates the build-up of the occupation number in absence of 
parametric amplification. 

6.1 No forbidden bands: formal proof 

The first order equations for the fermions ( |3.17| ) can be written in the following matrix form 



d , f / N 

*17 ~ M{t)a3 

dt 



fit) = (6.2) 



f{t)=[ ). (6.3) 



with cri_3 the Pauli matrices and 

It is straightforward to check that 

f\t)g{t) - fl{t)gi{t) + f;[t)g2{t) = constant (6.4) 



for an arbitrary pair of solutions f{t) and g{t) of the equation ( p.2[ ). This is a result of the conservation of 
probability in the Dirac equation, or equivalently the normalization condition on the Dirac spinors. Equation 



(3.2) has two linearly independent spinor solutions hi{t) , h2{t) with initial conditions: 



hm^i^ I ) ; hM={^ J )■ (6.5) 

Both solutions therefore obey 

h\{t)hi{t) = 1 , h\{t)h2{t) = 1 (6.6) 
Furthermore it follows from the above identities that 

h\{t)h2{t) ^ 0, (6.7) 

i.e, the linearly independent solutions are orthogonal at all times. 

Since Mlt) is periodic with period T, i.e, M{t + T) = M{t) then hi{t + T) and h2it + T) wiU also be solu- 
tions of ( |6.2[) , but then they must be linear combinations of the linearly independent solutions hi{t) , h2{t):, 



hi{t + T) = ahi{t) + (3h2{t) 
h2{t + T) = -fhi{t) + dh2{t) 

where a,(3,j,S are complex coefficients. The matrix of coefficients 



Z ? ) (6.8) 



7 5 

is called the monodromy matrix and is the important concept in Floquet theory. It represents an operator 
that evolves the solution in time by one period. 
The conditions 

h\{t + T)hi{t + T) = 1 ■ hl{t + T)h2{t + T) = 1 ; hl{t + T)h2{t + T) = (6.9) 
obtained above, lead to the following conditions on the coefficients 

lap + = 1 , |7|2 + |5|2 a*j + f3*S = 0, (6.10) 
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which in turn imphes that the monodromy matrix is unitary: 

M''M = 1. (6.11) 

This unitarity property of the monodromy matrix can be traced back to the conservation of probabihty 
of the solutions of the Dirac equation (since it is first order in time) as explicitly determined by the condition 

O- 

Floquet solutions are the eigenvectors of the monodromy matrix and the (logarithm) of the eigenvalues 
are the Floquet exponents. Floquet solutions therefore satisfy 

^±(i + r) =e'^±F±(t), (6.12) 

where in general, the exponents w± are complex. Their imaginary parts determine the growth (or decay) rate 
of the solutions and are responsible for parametric amplification. However, the unitarity of the monodromy 
matrix implies that the Floquet indices w± are real, hence Floquet solutions develop a phase upon time 
evolution during a period but their magnitude is constant. This then precludes forbidden bands and hence 
there is no parametric amplification of fermionic modes. Pauli blocking is a direct result of the Dirac fields 
obeying a first order evolution equation in time. This is also at the heart of the conservation of probability 
(or normalization of the Dirac spinors), which in turn determines that the monodromy matrix is unitary, 
hence Floquet indices are real. 

6.2 Perturbative argument: 

While the formal proof above unambiguously clarifies that there is no parametric amplification of fermionic 
modes with an oscillatory time dependent mass, we offer also a perturbative argument. We do this both to 
highlight the main result of the exact proof above, but also to illuminate why even when the Floquet indices 
are real and there are no forbidden bands, time evolution with phases leads to a build-up of the occupation 
number for some (resonant) values of the wave vectors. 

Consider the second order evolution equation for fi;p{t) (after rescaling by the scale M) 



dt^ 



+ + M^{t)+iM{t) 



M{t) = xo cos(m) 



(6.13) 



A perturbative solution is obtained by considering xo = £ < 1 a-nd writing the formal expansion fi^t) — 
fl^l(t) + ef^^l{t) + e^/-[^2(f) + • • • which leads to the following hierarchy of equations 



The zeroth order solution is of the form 



~^Mf['^.it) 



M'{t)f^'^Ht) 



B{p)e 



-ipt 



(6.14) 



3.15) 



with the coefficients A{p),B(p) determined by the initial conditions. The higher order inhomogeneous 
equations are solved in terms of the retarded Green's function 
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The general solution is given by 



G{p, t-t')^- sm\p{t ~ t')]Q{t ~ t'). (6.16) 
P 



fi'Ut)^ dt'G{p,t~t')l(^Ht') (6.17) 
Jo 

{i), 



with I*^*-' the inhomogeneity in the equation for the i-th order contribution. We note that ^(0) 



0; /{'^(O) = 0, therefore the initial conditions determine the constants A(p); B{p) to all orders. 

The first order solution f^p{t) exhibits a secular term, which grows linearly in time, when Q, = 2p. This 
obviously corresponds to the production of a fermion pair. However this secular term is purely imaginary: 
fi^pi^) "i-tXo/p- It corresponds to a renormalization of the phase of the order solution. The calculation to 
second order is lengthy but straightforward, with the following remarkable result: purely imaginary secular 
terms appear once again which grow in time, leading to a further contribution to the phases. The real part 
of the secular terms cancels out between the contribution (which involves the zeroth order solution ) and 
the contribution from iM which appears squared thus always (in both cases for /i,2) with opposite sign as 
the M^{t) contribution. Therefore, at least to second order, we find that the secular terms generated from 
resonances when f2 ~ 2p are purely imaginary, while the real part of the secular terms cancel exactly. Thus 
there is no exponential growth of the solution which would emerge from real secular terms in the perturbative 
expansion. We have not attempted a higher order calculation, but in view of the exact proof offered above, 
it is clear that the result holds to all orders. The lowest order perturbative solution also shows that /i and 



/2 develop opposite phases so that the occupation number given by eqn. (4.44) becomes of the form 



nq{T) cx sin('!i7r) (6.18) 

with Tu the real Floquet index. This explains the behavior of the occupation number displayed in the figures 
above, the oscillations and the saturation. Hence unlike the bosonic case in which parametric amplification 
results in an exponential growth of bosonic fluctuations, Pauli blocking at the level of the mode functions is 
a consequence of real Floquet indices that result in a bounded growth of the occupation number. 



7 Conclusions, conjectures, implications and further questions 

We have used an expansion in 1 /TV, where N is the number of fermion fields coupled to a scalar to understand 
the non-perturbative dynamics of this system. In this limit, the scalar fiuctuations are suppressed relative 
to the fermionic ones and we can consider mean field dynamics in a consistent approximation that, at least 
in principle, can be improved on. 

In the large N limit, the theory exhibits dimensional transmutation and symmetry breaking via the 
Coleman- Weinberg mechanism with an effective potential that is unbounded below and features a metastable 
minimum at the origin and two symmetric maxima. The bulk of our work consisted in obtaining the equation 
of motion for the expectation value of the scalar field, or mean field and comparing the dynamics due to the 
full equations of motion, including the backreaction of the fermionic fluctuations to that obtained from a 
renormalization group improved effective potential. 

The equation of motion obtained from the RG effective potential alone contains a Landau pole at a 
non-perturbative scale El oc e^/^~^ where g — y"^ /27r^ and y is the Yukawa coupling. We should note that 
since we are working within the large N expansion, this Landau pole cannot be considered to be an artifact 
of the perturbative expansion. The dynamics obtained from the effective potential alone would yield singular 
behavior when the amplitude of the mean field reaches the Landau pole. 

The renormalization of the full equations of motion leads to a running coupling constant that depends 
on time through the dependence on the amplitude of the mean field. This phenomenon is akin to the 
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dynamical renormalization found in reference ||2^. Furthermore we show that potential initial singularities 
are self-consistently removed and that the dynamics is smooth and free of Landau pole singularities, or 
discontinuities. 

The fully renormalized equations of motion including the backreaction of the fermionic fluctuations has 
two different regimes depending on whether the ultraviolet cutoff A <C El or A ^ E^. In the former case, if 
the initial value of the mean field is between the origin and the maxima, the mean field undergoes damped 
oscillations transferring energy to the fermionic fiuctuations. This energy transfer results in fermion pair 
production within a band of wavevectors determined by the mass of the scalar and the initial amplitude of the 
mean field. It is found that the backreaction from the fermionic fiuctuations, encoded in the renormalized 
fermion chiral condensate introduce a slowly varying renormalization of the parameters of the effective 
potential. 

When the absolute value of the initial value of the mean field is larger than the maximum of the effective 
potential, both the mean field and the fermion chiral condensate run away to the cutoff scale, at which point 
the evolution must be stopped. 

While it is commonly acknowledged that a theory that features a Landau pole only makes physical sense 
when the cutoff is below the scale of the Landau pole, we decided to study the dynamics even in the case in 
which A 3> El to find out if and how there are any singularities in the non-equilibrium dynamics associated 
with the Landau pole. Typical statements on the limit of validity of a theory with the Landau pole are 
mostly based on perturbative unitarity in S-matrix elements. 

Our study is definitely non-perturbative as we study the evolution for large amplitude mean field configu- 
rations. Whether perturbative statements on S-matrix elements provide a limitation to the non-perturbative 
approach is, to our knowledge, an open question. We therefore proceed to study the theory in this regime 
and to analyze the consequences without bias. 

In this case we find novel, remarkable behavior. When the amplitude of the mean field reaches the 
Landau pole there is a cancellation between the contribution from the effective potential and that from 
the renormalized chiral condensate that prevents a singularity in the dynamics at the Landau pole from 
developing. The ensuing evolution is continuous. As the amplitude of the mean field becomes larger than 
the Landau pole, it begins to oscillate with an amplitude that grows. We find that the time evolution of the 
mean field is well described by an effective potential that is quartic and upright. 

We have seen that the fermion occupation number almost saturates up to a wave vector of the order of 
the amplitude of the mean field. The resulting state is a very dense medium or plasm which can be described 
as a cold, degenerate fermi gas with a chemical potential of the order of the energy scale of the mean field. 

We have also provided an exact proof of the fact that the Floquet indices for the fermionic mode functions 
are real, therefore preventing unbounded parametric amplification. This is the manifestation of Pauli blocking 
at the level of the mode functions. A perturbative proof of this result served to illuminate the workings of 
the exact proof. 

Conjectures: We conjecture, at this stage, that the phenomenon of smooth dynamics beyond the 
Landau pole is somewhat related to a recent observation that in a medium the actual position of the Landau 
pole is actually shifted to much larger scales |Q. 

While we have studied the theory in a regime in which perturbative unitarity of S-matrix elements would 
suggest that the theory breaks down, we have done so non-perturhatively at the mean-field level. It is 
conceivable, that while the dynamics is smooth and there are no signals of pathologies, incompatibilities 
may lurk in some (non-perturbative) matrix elements or physical quantities that would cast doubt on the 
validity of the theory. To assess this possibility requires a deeper study of correlation functions beyond mean 
field and perhaps to next to leading order in the large N . This is beyond the goal of this article and an issue 
that deserves further scrutiny in its own right. 

There has previously been a related conjecture that Landau poles in some theories only appear in per- 
turbative calculations or in calculations of unobservablesp^ . 

Implications: There are many important cosmological implications of our results. While previous 
studies of reheating and preheating of fermionic theories have argued that fermions are not very efficient for 
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reheating because of Pauli blocking, in this work we have seen that in the regime when the cutoff is much 
larger than the Landau pole there is the possibility of abundant particle production. After the mean field 
crosses the Landau pole, the fermionic occupation numbers basically saturate up to a momentum scale of 
the order of the Landau pole. This in turn describes a dense medium with a large chemical potential of the 
order of the Landau pole. 

This could potentially be relevant in an extension of the Standard Model with top- like quark sector with 
a Yukawa coupling of order one. Since this theory will be part of an even larger theory (perhaps Grand 
Unified) the cutoff should be taken to be of the order the GUT scale. On the other hand with a vacuum 
expectation value of few hundred GeV's and a Yukawa coupling of order one, the position of the Landau 
pole is at a scale of a few TeV's. Thus the situation A ^ El would naturally arise in these type of scenarios. 

Another potentially interesting implication in cosmology would be the following: consider that after a 
phase transition the system is trapped in the metastable vacuum. A bubble type configuration, characterized 
by a collective coordinate[^ |2^, ^ will tunnel underneath the barrier of the effective potential|^ and exit 
on the other side of the maximum. The collective coordinate will then begin to roll down the potential hill 
until it reaches the Landau pole at which point the novel dynamics studied in this article will lead to fermion 
pair production and the production of a dense medium with a large chemical potential. We believe that 
these potential implications deserve further study and we expect to address these in future work. 

Further questions: An important aspect that we have not investigated is that of the new time scales 
that emerge from the dynamics. For example when the mean field oscillates with small amplitude in the first 
cases analyzed, there is a damping associated with the energy exchange with the fermion fluctuations. The 
time scale for damping is much larger than the typical oscillation frequency, despite the fact that the coupling 
is fairly strong. Another time scale is that associated with the growth of the amplitude of oscillations after 
the mean field crosses the Landau pole, until the amplitude saturates. 

We expect to apply the methods developed in||2^, ^to study these questions and relegate this investi- 
gation to future work. 
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